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Разработан двухэтапный метод повышения качества цветных изображений. Разработанный метод позволяет понижать уровень 
шумов и одновременно оптимизировать яркость и контраст изображения. Для удаления шумов применяется вейвлетная филь т- 
рация, для настройки яркости и контраста - генетический алгоритм. 


Введение 

Создание цифровых изображений с хорошей яр- 
костью, контрастностью и высоким уровнем дета- 
лизации является необходимым условием в различ- 
ных областях науки и техники. Среди них можно 
выделить: техническое зрение, дистанционное зон- 
дирование, биомедицинский анализ изображения, 
обнаружение неисправностей. Создание естествен- 
ных изображений или преобразование имеющихся, 
например, для улучшения восприятия информации 
в наглядной форме, является одной из важнейших 
задач обработки изображений. Методы, которые 
осуществляют такие преобразования, называют ме- 
тодами повышения качества изображения [1]. 

Большинство методов повышения качества изо- 
бражения, существующих на данный момент време- 
ни - это эмпирические или эвристические методы, 
зависящие от специфики изображения. Почти все 
эти методы требуют пользовательского взаимодей- 
ствия для настройки специальных параметров, что 
затрудняет их применение в автоматическом режи- 
ме. Наконец, существующие методы повышения 
качества изображения, как правило, подвергают 
преобразованию все изображение, в то время как во 
многих случаях необходимо преобразовать лишь ка- 
кую-то определенную область изображения. 

Описанный в [ 1 ] метод позволяет повышать ка- 
чество изображений автоматически, однако не 
всегда дает хорошие результаты. Кроме того, он 
применим только к полутоновым изображениям и 
требует достаточно много времени на обработку. 

В данной работе предпринята попытка разра- 
ботки полностью автоматического метода повыше- 
ния качества цветных изображений. Данный метод 
состоит из двух основных стадий, позволяющих 
исправить основные недостатки цифровых изобра- 
жений: удаление шумов и оптимизация яркости и 
контраста. 

1. Коррекция яркости и контраста 

1.1. Генетический алгоритм 

Генетический алгоритм (ГА) представляет собой 
метод оптимизации, основанный на концепциях 
естественного отбора и генетики. В этом подходе 
переменные, характеризующие решение задачи, 


представлены в виде генов в хромосоме [2]. ГА, 
оперируя конечным множеством решений, генери- 
рует новые решения как различные комбинации 
частей решений популяции используя такие опера- 
торы, как отбор (селекция), рекомбинация (крос- 
синговер), мутация. Новые решения позициониру- 
ются в популяции в соответствии с их положением 
на поверхности исследуемой функции [2]. 

Базисом ГА являются: 

• особь с определенным набором хромосом; 

• начальная популяция особей; 

• набор операторов для получения новых реше- 
ний; 

• целевая функция для отбора наиболее приспо- 
собленных решений. 

Основными операторами генетического алго- 
ритма являются: 

• селекция (отбор наиболее приспособленных 
решений на основе целевой функции); 

• кроссовер (скрещивание особей для получения 
новых наборов хромосом); 

• мутация (случайное изменение одного или бо- 
лее генов особи). 

В разработанном алгоритме каждая особь имеет 
4 хромосомы, которые представлены беззнаковы- 
ми длинными целыми числами (32 бита). 

1.2. Ядро улучшения 

Методы локального улучшения основаны на 
применении функции, зависящей от распределе- 
ния яркости в окрестности каждой точки изобра- 
жения. Одним из примеров таких методов является 
адаптивное выравнивание гистограмм, показавшее 
хорошие результаты в обработке медицинских изо- 
бражений. Однако этот метод использует ядро, тре- 
бующее больших вычислительных затрат. 

Ядро улучшения - это преобразование Г, при- 
меняемое к каждой точке изображения с координа- 
тами (х, у). Входным параметром для этого преоб- 
разования является исходная яркость точки полу- 
тонового изображения Дх, у), выходное значение - 
итоговая яркость точки §(х, у). Пусть Ѵ ще и - со- 
ответственно вертикальный и горизонтальный раз- 
мер изображения. Тогда преобразование Т опреде- 
лим следующим образом [ 1 ] : 
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где §(х,у) и Ах, у) - полученная и исходная яркость 
точки с координатами (х,у), М - глобальное сред- 
нее значение яркости изображения, т(х,у) - сред- 
нее значение яркости в окрестности пхп, а(х,у) - 
значение среднеквадратического отклонения в 
этой же окрестности, а, Ъ, с и к - параметры ядра 
улучшения, наиболее подходя щ ие значения кото- 
рых и определяются при помощи генетического ал- 
горитма. 

В ходе проведения исследований было выявле- 
но, что применение среднеквадратического откло- 
нения в чистом виде приводит к не очень хорошим 
результатам, поскольку в однотонных областях его 
значение близко к нулю, что приводит к сильному 
осветлению. В итоговом алгоритме при вычисле- 
нии среднеквадратического отклонения в окрест- 
ности пхп к ней добавлялась черная или белая 
область такого же размера. Добавление черных 
областей позволяет добиться равномерного освет- 
ления слишком темных изображений, добавление 
белой области - наоборот, позволяет затемнить 
слишком светлые. 


1.3. Критерий улучшения качества изображения 

Для автоматизации процесса улучшения каче- 
ства изображения необходимо определить объек- 
тивный критерий для оценки качества улучшения 
изображения. Этот критерий используется для по- 
строения целевой функции генетического алгорит- 
ма и включает в данном алгоритме 3 составляющих: 

1 . Количество и интенсивность краевых пикселей. 

Улучшенное изображение имеет большее число 
краевых пикселей с более высокими значениями 
интенсивностей по сравнению с исходным изобра- 
жением. Для того, чтобы выявить эти характери- 
стики, применяется так называемый детектор края 
изображения. В данной работе был применен опе- 
ратор Собеля [5]. 

Оператор Собеля [5] производит измерение 
двумерного пространственного градиента на изо- 
бражении и выявляет области с большим значени- 
ем этого параметра. Эти области и соответствуют 
краям. Как правило, он используется для оценки 
модуля градиента в каждой точке черно-белого 
изображения. 

Работа этого оператора основана на наложении 
на каждую точку изображения двух масок враще- 
ния. Эти маски представляют собой две ортого- 
нальные матрицы размерностью 3x3 (рис. 1). Эти 
маски выявляют границы, расположенные верти- 
кально и горизонтально на изображении. При раз- 


дельном наложении этих масок на входное изобра- 
жение можно получить оценку градиента по каж- 
дому из направлений (обозначим их О, и О у ). Ко- 
нечное значение градиента получается по формуле 

С(х,У) = ^ х 2 ( х ’У) + С у 2 ( х ’У'>- 
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Рис. 1. Маски оператора Собеля 


При помощи данного оператора подсчитывает- 
ся число краевых пикселей г](І) и их суммарная ин- 
тенсивность Е(І): 

Е(І) = ^С(х,уУ, 

X у 

гі(і) = ^н( х ,уУ, 

х,у 


рі(х,у) 


Гі ,С(х,у)>а 
[0,С(х,_у) < а 


здесь а - некое предельное значение. 

2. Мера энтропии изображения. 

Необходимо оценить, каково количество пик- 
селей с различным уровнем яркости присутствует 
на изображении. 

Для этого необходимо построить гистограмму 
изображения, количество уровней которой соответ- 
ствует количеству градаций яркости. Предельные 
значения уровня с индексом і записываются как А, и 
Д для /е{0,...,255}. На основе гистограммы вводит- 
ся мера энтропии улучшения изображения I, как 

Я(/) = 2>( 0; 


ѵ(і) 


Г-Ѵ.-ІО е 2 (ѵ,), V,. 
[О, ѵ і = 0 


(*) 


где у, - частота пикселей, имеющая уровни града- 
ций яркости между граничными значениями А, и Д. 
Отметим, что мера в (*) это не реальная энтропия, 
т. к. ѵ, - скорее относительные частоты, чем веро- 
ятности. 

3 . Уровень адаптации к зрению человека по яркости. 

Для человеческого зрения наиболее приятны 
изображение, среднее значение яркости на кото- 
рых равно середине всего диапазона возможных 
яркостей [3]. Поэтому в целевую функцию был вве- 
ден следующий множитель: 

7-7 /7 

10 = 1-0,3 

7/2 

шах 

Здесь Х шах и 7 - максимальное и среднее значе- 
ние яркости на изображении. 

Таким образом, целевая функция должна быть 
пропорциональна Е(І), Н(Г) и 10. Была вы- 
брана функция, имеющая следующий вид: 
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Р{х) = 1п(ВД/)) + е) 77 ( 2 е" (І ЧѲ, 

Н . УѴ 

5іхе зіге 

где Н !і№ и ]Ѵ !І№ - соответственно вертикальный и го- 
ризонтальный размер изображения, е=2,718. 

Данная функция очень близка по своему виду к 
функции, использованной в [1]. 

1.4. Улучшение цветных изображений 

Рассмотренный выше метод позволяет улучшать 
изображения с одной составляющей яркости - по- 
лутоновые изображения. Для повышения качества 
изображений, имеющих три цветовые составляю- 
щих яркости пикселя, они первоначально преобра- 
зуются в полутоновые. 

После этого все преобразования производятся 
над черно-белым изображением. 

По окончанию работы алгоритма изображение 
снова переводится в цветное путем умножения 
каждой компоненты яркости исходного изображе- 
ния на соответствующий коэффициент. 

2. Удаление шумов 

2.1. Вейвлет-преобразование 

Непрерывное вейвлет-преобразование (Сопіі- 
П 110118 \Ѵаѵе1ег Тгашіогт, С\ѴТ) описывается ура- 
внением (3), где х(і) - анализируемый сигнал, у/(1) 
- вейвлет-функции, получаемые из материнского 
вейвлета у/({) путем сдвига и масштабирования (ра- 
стяжения или сжатия) [4]. 

ѵМ 

Параметр сдвига г определяет положение вей- 
влет-функции относительно сигнала, т. е. он явля- 
ется временной характеристикой вейвлет-преобра- 
зования. Параметр масштабирования 5 несет ча- 
стотную информацию. Материнский вейвлет, ис- 
пользуемый для получения базисных функций, вы- 
бирается на основе каких-либо необходимых 
свойств этих функций. Масштабирование по вре- 
мени растягивает либо сжимает сигнал. Низкоча- 
стотные составляющие несут общую информацию 
о сигнале, высокочастотные - информацию о де- 
талях сигнала. 

Вейвлет-ряды получаются путем дискретиза- 
ции вейвлет-преобразования по времени, которая 
позволяет производить вейвлет-преобразование с 
использованием компьютеров. 


Однако, вейвлет-ряды это лишь выборочная 
версия непрерывного вейвлет-преобразования и их 
вычисление может требовать больших затрат време- 
ни и ресурсов. Дискретное вейвлет-преобразование 
(Оінсгсге \Уаѵе1еІ Тгашіогт, 0\ѴТ) разработано для 
быстрого вычисления вейвлет-преобразования. Оно 
просто в выполнении и позволяет уменьшить время 
вычисления и количество требуемых ресурсов. 

В С\ѴТ - п реобразо ва н и и сигналы анализируют- 
ся с использованием набора базисных функций, 
получающихся друг из друга путем сдвига и мас- 
штабирования. В случае ШѵТ-преобразования 
масштабное временное представление сигнала по- 
лучается путем цифровой фильтрации. Анализиру- 
емый сигнал проходит через фильтры, на которых 
происходит «обрезание» определенных частот с 
различными масштабами [4]. 

2.2. Удаление шумов при помощи ОѴѴТ 

Большинство фильтров, основанных на вейвлет- 
преобразованиях, имеют одинаковую схему (рис. 2). 

Рассмотрим все части данной схемы примени- 
тельно к разработанному алгоритму. 

Первоначально изображение, представленное в 
формате КС В, переводится в цветовое простран- 
ство УЦѴ. Компоненты И и V нс подвергаются из- 
менению во время фильтрации. Входным сигналом 
для системы является яркостная составляющая 
изображения (У) [7]. 

У = 0,2997? + 0,578С + 0,114Д 

V = 0,16787?-0,3313С + 0,5В; 

V = 0, 57? - 0, 4 1 87 С + 0, 08 1 3 # 

Вейвлет-преобразование в данном случае - это 
двумерное дискретное вейвлет-преобразование [7], 
основанное на материнском вейвлете «3-5», имею- 
щем 3 высокочастотных и 5 низкочастотных коэф- 
фициентов (рис. 3). 

Фильтрация в нашем случае представляет собой 
удаление из полученного набора вейвлет-коэффи- 
циентов высокочастотных компонент, которые 
включают и шумы. Пользователь сам может опре- 
делить порог при выполнении этой операции. 

После обратного вейвлет-преобразования изо- 
бражение переводится обратно в цветовое про- 
странство КС В по формулам: 

К = У + 1,3707 -(К -128); 

О = У -0,3365-(і/-128) -0,6982(К-128); 

В = У + 1,7324 •({/ -128). 


Входной 

сигнал 



Выходной 

сигнал 


Рис. 2. Обобщенная схема фильтрации 
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Рис. 3. Материнский вейвлет, использованный в разрабо- 
танном алгоритме 

Первоначально планировалось производить 
удаление шумов до коррекции яркости и контра- 
ста, поскольку во время этой стадии обработки 
присутствующие на изображении шумы усилива- 
ются. Однако, как показали исследования, данный 
вид фильтрации снижает контраст изображения, и 
при этом могут потеряться какие-то детали, слабо 
различимые на исходном изображении. Поэтому 
было решено производить удаление шумов после 
подстройки яркости и контраста. 

3. Результаты работы 

Итогом исследований явилось создание прило- 
жения, реализующего данный метод улучшения 
изображений. Это приложение обрабатывает фай- 
лы формата *.Ътр размером до 512x512 пикселей. 
Приложение было применено к набору тестовых 
изображений. Оказалось, что данный метод хоро- 
шо проявляет себя на двух классах изображений: 

1. Слишком светлые или темные изображения. 
Данный метод позволяет добиться осветления 
слишком темных участков и затемнения слиш- 
ком светлых. Кроме того, он позволяет про- 
явить объекты, которые на исходном изображе- 
нии были практически неразличимы (рис. 4, 5). 

2. Затуманенные изображения. Применение данно- 
го метода к затуманенным или нечетким изобра- 
жениям позволяет повысить их четкость, а также 
в некоторой мере снять затуманенность (рис. 6, 7). 
Реализованный генетический алгоритм при раз- 
мере популяции 20 особей находил достаточно хо- 
рошее решение за 5 поколений. При этом обработка 
изображения в среднем занимает 10.. .20 с (напри- 
мер, для изображения размером 498x326 пикселей 
обработка заняла 16 с на компьютере с процессором 
Іпіеі Репіішп 4 2,4 ГГц и ОЗУ объемом 512 Мб). 

Заключение 

Описанный метод позволяет бороться с основ- 
ными дефектами изображений - шумами, низким 
контрастом, неоптимальной яркостью и т. д. Вся 
обработка происходит в автоматическом режиме и 


занимает относительно немного времени по срав- 
нению с другими методами, использующими те же 
принципы. 

Работа поддержана грантом № 06-08-00840 Рос- 
сийского фонда фундаментальных исследований. 




Рис. 4. Сопоставление изображений: а) до обработки, Ь) после 



а Ь 

Рис. 5. Сопоставление изображений: а) до обработки, Ь) после 
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а Ь 

Рис. 6. Сопоставление изображений: а) до обработки, Ь) после 



Сопоставление изображений: а) до обработки, Ь) после 



с ё 

Рис. 8. Использование всех шагов алгоритма для коррекции изображений. Изображение: а) исходное, Ь) зашумленное, 
с) улучшенное, б) после удаления шумов 
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